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c3 ■ Abstract 

We investigate spectral functions extracted using the Maximum Entropy Method from cor- 
relators measured in lattice simulations of the (2+l)-dimensional four-fermion model. This 
model is particularly interesting because it has both a chirally broken phase with a rich spec- 
trum of mesonic bound states and a symmetric phase where there are only resonances. In 
the broken phase we study the elementary fermion, pion, sigma, and massive pseudoscalar 
meson; our results confirm the Goldstone nature of the vr and permit an estimate of the 
meson binding energy. We have, however, seen no signal of a ^ vrvr decay as the chiral 
limit is approached. In the symmetric phase we observe a resonance of non-zero width in 
qualitative agreement with analytic expectations; in addition the ultra-violet behaviour of 
the spectral functions is consistent with the large non-perturbative anomalous dimension for 
fermion composite operators expected in this model. 
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1 Introduction 



The Gross-Neveu model in d = 3 spacetime dimensions (GNM3) has been the object of much 
analytic and numerical study in recent years. Its Lagrangian density is 

~ iJi(^^ + m + -^^{cr + ij57i)^iji + ^{a^ + 7r^), (1.1) 

where the index i runs over Nf fermion flavors and in the second line we have introduced 
scalar a and pseudoscalar vr auxiliary boson fields. Apart from the obvious numerical ad- 
vantages of working with a relatively simple theory in a reduced dimensionality there are 
several features which make GNM3 interesting for the modelling of strong interactions 

• For sufficiently strong coupling > g1 it exhibits spontaneous chiral symmetry break- 
ing implying dynamical generation of a fermion mass Mj, the pion field vr being the 
associated Goldstone boson. A separation of scales -C Mf is possible. 

• The spectrum of excitations contains both "baryons" and "mesons", namely the ele- 
mentary fermions / and the composite // states; 

• For 2 < (i < 4 there is an interacting continuum limit at a critical value of the coupling, 
which for = 3 has a numerical value g1/a ~ 1.0 in the large- iVj limit if a lattice 
regularisation is employed There is a renormalisation group UV fixed point at 
g"^ = g1i signalled by the renormalisability of the 1/A^/ expansion ||l[, entirely analagous 
to the Wilson-Fisher fixed point in scalar field theory. 

• Numerical simulations with baryon chemical potential ;U 7^ show qualitatively correct 
behaviour, in that the onset of matter occurs for /x of the same order as the constituent 
quark scale Mf 0, rather than for /i ^ "^77/2, which happens in gauge theory simu- 
lations with a real measure det (M^M) because of the presence of a baryonic pion in 
the spectrum. This makes GNM3 an ideal arena in which to test strongly interacting 
thermodynamics 0. 

Let us briefiy review the physical content of the model as predicted by the large-A^j 
approach [|^, |^. For g'^ > g^ the fermion has a dynamically generated mass Mf given, up to 
corrections of order l/Nf, by 

Its inverse defines a correlation length which diverges as [g"^ — g1)~'^ with critical index 
z/ = 1 + 0{1/Nf). In addition as a result of // loop corrections the a and vr fields acquire 
non-trivial dynamics, the inverse a propagator being given as a function of d to leading order 
in 1/A^/ by 



1 (47r)i M^-^ 

?2r(2-f) (P + 4M2)F(1, 2 - f ; I; -^) 



= .....ZL (1-3) 
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Immediately we see the difference between this model and QCD. For k"^ <^ Mj F ^ 1, 
implying that to this order the a resembles a weakly-bound meson of mass = 2Mf; 
however, the hypergeometric function in the denominator indicates a strongly interacting 
// continuum immediately above the threshold 2Mf. This implies that if truly bound, its 
binding energy is 0{1/Nf) at best (to our knowledge there have so far been no analytic 
calculations), implying little if any separation between pole and threshold. Since all residual 
interactions are subleading in 1/Nf, we surmise that all other mesons are similarly weakly 
bound states of massive fermions, and hence effectively described by a two-dimensional 
"non-relativistic quark model" . A recent study of mesonic wavefunctions in GNM3 provides 
evidence for this picture 0. In an asymptotically-free but confining theory like QCD, by 
contrast, one expects isolated poles and/or resonances, corresponding to relativistic bound 
states in the channel in question, which are well separated from a threshold to nearly-free 
quark behaviour which sets in at typically 1.3 - 1.5 GeV 0. 

The exception to this rule is the pion. The Lagrangian ( |1 . 1| ) can be defined with either 
a continuous U(l) or discrete Z2 chiral symmetry, the latter case being realised by setting 
the 71 field to zero. In the case of U(l) chiral symmetry, for m = and > g1 the pion 
propagator D^^ is given by a similar expression to ( [L.3|) with the factor (fc^ -|- 4M|) replaced 
by fc^; the massless pole demonstrates that vr couples to a Goldstone mode. For m > 0, we 
expect by the usual PCAC arguments that the vr acquires a mass m,r oc y/m, and that the 
ratio rriTr/Mf can be tuned to be arbitrarily small. In particular, once it is less than unity 
the a becomes unstable with respect to decay into 2tt. Note, however, that the Goldstone 
mechanism in GNM3 is fundamentally different from that in QCD. In GNM3 the diagrams 
responsible for making the pion light are flavor- singlet chains of disconnected // bubbles 0. 
The non-singlet connected // diagram which interpolates the pion in QCD corresponds in 
GNM3 to a pseudoscalar state with mass 0(2M/). 

For g"^ < g1 the model is chirally symmetric, and hence all states are massless, as m ^ 0. 
In this limit and D^^ coincide, and in the large-A^j limit neither has a pole on the 
physical sheet |]l|. The auxiliary fields in this case do not interpolate to a stable particle. A 
dimensionful scale is still defined, however, by the width of a resonance in // scattering 
in these channels; this diverges as {g1 — g"^)^" with the same exponent u 0. 

It is clear that despite its simplicity GNM3 exhibits phenomena such as resonances, decays 
and multi-particle continua which are not easily analysed using the traditional techniques of 
single- and multi-exponential fitting to Euclidean correlators developed for quenched QCD. 
This was recognised in early studies, which attempted fits inspired by the large-A^j forms 
of Dfj in both chirally broken and symmetric phases, with ambiguous results 0. A more 
systematic approach, however, is to focus on the spectral density function p{u!), defined 
implicitly via the Euclidean timeslice meson correlator C{t) by 



Here J is a local fermion bilinear ipTip which in principle projects onto all physical states 
consistent with a given set of quantum numbers. All information about bound states, reso- 
nances and particle production thresholds as a function of energy u is contained in p. The 
procedures for fitting lattice-generated data to date have assumed Ansdtze for p such as one 




(1.4) 
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or more bound state poles of the form S{uj — M), or perhaps a free particle continuum above 
some threshold 0. However, more recent works have attempted ab initio calculations of 
p(co') 0. This is a difficult problem: the inversion of ( p..4| ) is ill-posed since p{uj) is a 



continuous function whereas lattice simulations only yield C{t) for a discrete, finite set of 
points, and moreover with some statistical uncertainty. The approach adopted in Refs. |0 
is to apply the Maximum Entropy Method (MEM) which attempts to fit p{uj) subject to 
reasonable assumptions of smoothness and stability with respect to small variations in the 
input data. 

In this paper we present results from a study of spectral functions extracted from numer- 
ical simulations of GNM3 using MEM techniques. To our knowledge this is the first such 
study beyond the quenched approximation. Our goal is to explore some of the features de- 
scribed above which distinguish GNM3 from quenched QCD. In this regard it is worth noting 
that because the two most important mesonic channels, a and vr, are represented by bosonic 
auxiliary fields, the correlation functions in these channels automatically include the discon- 
nected diagrams which are so expensive to calculate in QCD; in GNM3, by contrast, these 
can be measured with high statistics relatively cheaply. We will also examine the fermion 
and non-singlet pseudoscalar (PS) channels. As surveyed above, simulations of GNM3 offer 
the freedom to vary the phase of the theory (by varying sgn(5f^ ~ 9l))i correlation length 
(by varying \g'^ — (yf^|), the symmetry group (by including or omitting vr), the ratio niT^/Mf 
(by varying m), and the interaction strength (by varying Nf) - in the current study we will 
exploit most of these opportunities. In future work we plan also to study the model with 
both non-zero temperature T and baryon chemical potential p. 

In Sec. ^ we survey MEM and explain our implementation of it. Sec. |^ outlines some 
theoretical expectations related to p{(jj) in GNM3 based on the large- A^/ approach, and Sec. ^ 
details the lattice formulation and numerical simulations. Our results are presented in Sec. ||, 
and conclusions in Sec. |6|. 



2 The Maximum Entropy Method 

The theoretical basis for MEM is Bayes' theorem in probability theory ||TT||: 

where P[X|y] denotes the conditional probability of X given Y . In terms of the lattice data 
D, spectral function p and all a priori knowledge if, Bayes' theorem reads: 

lA^I J P[D\H] ^ ' 

P\P\pH\ is known as the likelihood junction and is the equivalent of the familiar in the 
least squares method [jl2l . For a large number of Monte Carlo measurements of a correlation 
function, the data D are expected to obey a gaussian distribution according to the central 
limit theorem, 

P\D\pE\ = ^e-^^P\ (2.3) 
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L[p] 



Nt 



(2.4) 



where the normahsation factor = (27r)^ Vdet C and Nt is the number of temporal points. 
Lattice data averaged over A^^fg gauge configurations D(t), the covariance matrix Cij, and 
the propagator constructed from the spectral function p using the lattice kernel K{uj, t) are 
defined by 



cfg 



N, 



cfg m=l 
1 



iVcfg(A^cfg 

K(t, uj)p{ijj)duj 



-Y.{D"^{U)-D{U)){D-{t,)-D{t,)l 

m=l 



(2.5) 

(2.6) 
(2.7) 



In all our work we use a lattice kernel defined as exp(— tut). 

The factor P[p\H] appearing in the numerator of ( |2.2| ) is the prior probability, which is 
written in terms of the Shannon- Jaynes entropy S[p] |jl3| for a given default model po{uj). The 
default model is usually chosen to be the spectral function for a non-interacting two-particle 
continuum; for meson states we have Pq{uj) oc o;"^"^ (see Sec. ||). The final result, however, 
should be insensitive to the choice of po- The entropy S[p\ < and becomes zero only when 
p{uj) = poiuj): 



P[p\Hapo] 
S[p] 



Zs{a) 



^aS[p] 



E 



p{uj) - poiuj) - p{uj) In 
Pe - Poe - Pe In 



.Poi, 



Poiuj) 



duj 



(2.8) 
(2.9) 
(2.10) 



where Eq. (|2.10|) results from discretising the cj-axis into N^^ bins of width Acu, and the 
normalisation factor Zs = {^)^ ■ Note that two extra parameters previously implicit in 
H have been written in explicitly; a is a real positive parameter and po{uj) a real positive 
function. This form of entropy leads to a positive semi-definite spectral function in MEM. 
In our work we use A^^^^ = 600 and Acua = 0.01. 
Combining Eqs. ( p.3D and (|2.8| ) gives 



P[p\DHapo\ oc 
Q ^ 



1 



oQ[p] 



ZsZl 
aS — L. 



and the condition satisfied by the most probable spectral function Paiyj) is 

6Q 



Sp{uj) 







(2.11) 
(2.12) 

(2.13) 



P=Pc, 
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The parameter a is in effect a relative weighting between the entropy S and the hkehhood L, 
and there are three different ways to deal with it. The value a = a can be chosen which either 
gives = '^i or maximises P\a\DH p^] these methods are known as classic and historic 
||13| respectively. Alternatively, a weighted average over P[a\DHpo] can be performed; this 
is known as Bryan's method |M and is the one we adopt: 



dap^{uj)P[a\DHpo], (2.14) 

where a^m and amax are chosen to satisfy 

P[a^m,max\DHpo] = 0.01P[a\DHpo]. (2.15) 



2.1 Testing MEM 



To test our implementation of MEM, we studied an idealised QCD spectral function in the 
charged p- meson channel P, H: 



2 

TT 



1 



where the pole residue Fp = fpurip is defined by 

{Q\d-ipu\p)tp = y2/pmje^. 



(2.16) 



(2.17) 



being the polarisation vector. The following energy-dependent width is chosen with a 
^-function included to give the correct threshold behaviour of the p — vrvr decay 



rp(^) 



^mJl 
487r ^ I 



4m: 



2\ 3/2 



e{uj - 2m^). 



(2.18) 



The values of the parameters input into ( p.l6| ) are taken to be 



rrir 



9p-n 



0.77, 

= 5.45, /„ 



0.14, ujQ 
0.184, 5 



1.3, 
0.2, 



a. 



0.3. 



(2.19) 



where the numerical values of the first three parameters are in GeV. 

Test lattice correlator data were constructed from the spectral function using ( |2.7| ) . Gaus- 
sian noise with variance o"(tj) = 6Dj„(tj)ti was added to this data to simulate the effect of 
decreasing signal-to- noise ratio with temporal separation For simplicity we use a diagonal 
covariance matrix, which thus neglects correlations between different t. The default model 
used is po{uj) = mouj"^, motivated by the asymptotic behaviour of The parameter mo is 



chosen to be lim,, 



0.0277. We set cj. 



6GeV, Aco = lOMeV and N,., 



600, 



and vary the noise parameter b from 0.1 to 0.001. Fig. [l| shows a comparison between pj„ 
and Pout for various b. As expected, decreasing b leads to a better agreement between input 
and output spectral functions. 



6 




Figure 1: Comparison of pj„ and pout for idealised data in the p-channel with Nt = 32. 

3 Theoretical Preliminaries 

Our main focus will be the mesonic Euclidean timeslice correlation functions defined in 
Eqn. ( |1.4| ). With this definition, if J couples to a stable (ie. zero width) bound state of 
mass M with strength A (ie. {0\J\k,M) = A), then p{uj) = {\A\'^ /2M)6{uj - M). Since in d 
spacetime dimensions the engineering dimension [J] = d — 1 and M)] = 1 — |, it is readily 
checked that the combination p(uS) j uj'^^'^ is dimensionless. This also motivates the use of the 
default model p^iyj) oc cj'^"^, which corresponds in configuration space to the propagation of 
free massless fermions; ie. C{t) oc t~^'^~^\ For an asymptotically-free theory such as QCD we 
expect lim^^ooP(i^) = Po{^), as illustrated in Fig. ^ however since GNMs's UV behaviour 
is described by a renormalisation group fixed point with non-vanishing interaction strength 
m, H this is not a constraint in the current study. 

The asymptotic form of p{uj) is easiest to analyse in the symmetric phase < of the 
model, where we have a large- A^/ prediction 0. In the scalar channel, the momentum 
space propagator 

where 2 < d < A and /i is a dimensionful scale which increases as {g"^ — f?^)", ie. as an 



inverse correlation length. For d = 3 this implies []16 



C^{t) oc p / dk- = p / diu— -e-'^\ (3.2 



and hence the large- A^/ prediction 



k + p Jo uj'^ + p- 



p„{u) oc (3.3) 
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In the asymptotic regime we thus have p„ puyioj) oc uj~^ rather than poioj) oc u. This is 
a consequence of the large non-perturbative anomalous dimension 77^^ = d — 2 acquired by 
the scalar density at the UV fixed point 0, which relates the asymptotic forms via 

Puvi'^) Pq{uo)uo^'^'^^^ . (3.4) 

At smaller energy scales we interpret p as describing a resonance whose central position and 
width are both 0{p) and hence increase as the coupling g'^ is reduced. A second prediction 
of ( p.3|) is that the dimensionless combination p{uj)/uj tends to a constant in the IR limit 
0. 

Another situation of interest is the possibility of a decay in the chirally broken phase. 
Denote the physical fermion mass by Mj; the a is then expected to be a weakly bound state 
of mass Mo- ^ 2Mj whereas, for the case of a continuous chiral symmetry, the pion mass 
m^r may be much smaller. If 2m,r < Mo-, the decay a 27r is allowed and should show 
up as a threshold in the scalar spectral function. This should be a good warm-up exercise 
for studying the physical decay p ^ 27r in QCD; as well as the computational saving, an 
important technical consideration in the present case is that unlike in QCD the two pions 
can be produced in a state of zero relative momentum. 

Let us first derive an expectation for the form of the threshold using the 1/iV/ expansion. 
The contribution of the two pion intermediate state to the a correlator is shown diagra- 
matically in Fig. |^. To leading order in 1/Nf, using the conventions of Sec. 2 of the a 




Figure 2: Contribution to D^j from a 27r intermediate state. 



propagator is given by Eqn. (|1.3| ) where for momenta k <C Mf the hypergeometric func- 
tion in the denominator may be approximated hj F ^ 1. We will assume that for bare 
fermion mass m > 0, the pion propagator is given by the same expression with the factor 
{k^ + 4M|) in the denominator replaced by (fc^ -|- m^). The vertex Fo-TTTr is assumed to arise 
from a single fermion loop as indicated in Fig. ^ It is identically zero if chiral symmetry is 
unbroken. Using the bare vertex —g/yjN^, it is straightforward to show 

V ~ —R ^ {'\ ^\ 

CTTTTT '--^ (TTTTT \^*^/ 

where Gotttt is a dimensionless d-dependent constant. 
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With these components in place it is now possible to calculate D^j including the effects 
of the two pion intermediate state. Specialising to d = 3, we find 



fc2 + 4M2- i-tan- , 



(3.6) 



Besides the pole at k"^ ~ —AMj, there is now a contribution at 0{1/Nf) to the timeslice 
correlation function given by 



Gl^K-n rdk e'« 



The two pion threshold manifests itself via branch cuts in the inverse tangent running from 



k"^ = —Ami out to ±ioo. Approximating k"^ <^ Mf as before we integrate around the cut in 



from which we read off 



the upper half plane to obtain 

Eqn. ( p.9[ ) predicts that as well as a pole at a; ~ 2M/, there should also be a spectral 
feature at a; = 2mT, whose strength scales as {NfAdfrn-,,)'^] this is in principle testable by 
varying the simulation parameters Nf, g"^ and m. On a finite volume it will, however, prove 
difficult to study the detailed form of the spectral function above threshold. This is because 
the number of modes into which the a can decay is strictly delimited by the allowed pion 



wavevectors k^^ = 27rn/Ls, where n has integer-valued components, and 2ym| + k"^ < M„. 
The optical theorem, however, implies that the only intermediate states which can contribute 
to p{uj) are possible decay modes of the a; we infer that on a finite lattice, the uj~^ shape 
predicted by ( |3.9| ) is replaced by a set of 5-functions, each arising from an allowed /c^. With 
imperfect (ie. finite) statistical data, however, it is possible that under MEM these isolated 
poles will blend into a continuum of approximately the correct shape. 



4 Simulations 

The fermionic part of the lattice action we have used for the semi-bosonized GNM3 with 
U{1) chiral symmetry is given by [Q 

= E \ T.Xi{x)M.yXiiy) + lT.X^{x)x^{x)[J2a{x)+^e{x) > (4-1) 

i=l \x,y ° X (x,a;) {x,x) J 

where Xi ^^id Xi are Grassmann-valued staggered fermion fields defined on the lattice sites, 
the auxiliary fields a and vr are defined on the dual lattice sites, and the symbol {x, x) denotes 
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the set of 8 dual lattice sites x surrounding the direct lattice site x. The fermion kinetic 
operator A4 is given by 

■^x,y = T^Y.V'^i^) [^y,^+i> - ^y,x-i>] + ^^x,y, (4.2) 

where ri^{x) are the Kawamoto-Smit phases (— 1)^°"' '"^"-i^ and the symbol e(x) denotes 
the alternating phase i)^o+^i+^'2. The auxiliary fields a and vr are weighted in the path 
integral by an additional factor corresponding to 

Saux = ^,T.[As^) + As^)]- (4.3) 

The simulations were performed using a standard hybrid Monte Carlo (HMC) algorithm 
without even/odd partioning, implying that simulation of N staggered fermions describes 
Nf = AN continuum species the full symmetry of the lattice model in the continuum 
limit, however, is U{Nf/2)v ® U{Nf/2)v ® U{1) rather than U{Nf)v ® U{1). At non-zero 
lattice spacing the symmetry group is smaller still: U {N f / 4:)v ®U {Nf / 4:)v ®U {!) . In the Z2- 
symmetric model the vr fields are switched off and M becomes real so that real pseudofermion 
fields can be used. In this case staggered fermions describe Nf = 2N continuum species. 
Further details of the algorithm and the optimisation of its performance can be found in 

Bi- 

Using point sources we calculated the zero momentum fermion (/) correlator at different 
values of the coupling (3 = 1/ g"^. In order to compare MEM to conventional spectroscopy we 
also estimated the fermion mass using a simple pole fit using the function 

Cfit) = Af[e-^''i' - (-l)*e-^^/(^*-*)]. (4.4) 

Similarly, the zero momentum auxiliary vr correlator was measured and its mass estimated 
using a cosh fit. The mesonic correlators are given by: 

CM{t)= J2 <f(xi)<I>(x2)W^M(x)G(x,t;xi,0)Gl-(x,t;x2,0), (4.5) 

X,Xl,X2 

where G is the lattice fermion propagator and ^^/(x) a phase factor which picks out a 
channel with particular symmetry properties i.e. Wm{'x) = e{x) for the S channel and 
Wm(x) = 1 for the PS channel. The function $(x) is either a point source 5x,(o,o) or a 

staggered fermion wall source Ylm,n=o ^:x.,(2m,2n) lOl- simulations we used point 

sinks. These correlators were fitted to a function CM{t) given by 

CMit) = A[e-^'''' + e-*^«(^*-*)] + i(-l)*[e-*^^^* + e-^''''^^'-% (4.6) 

Note that composite operators made from staggered fermion fields project onto more than 
one set of continuum quantum numbers. The first square bracket represents the "direct" 
signal with mass Mm and the second an "alternating" signal with mass Mm- Continuum 
quantum numbers for various mesonic channels are given in ^ - in this study we focus on 
the PSdirect channel, with = 0~ . Although expected to be the tightest bound meson since 
it is the only one for which s-wave binding is available, as stressed in 0, ^ this state does 
not project onto the Goldstone mode in the broken phase. 
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5 Results 



5.1 The TT, / and PS Channels in the Broken Phase 

We first discuss results from tlie cliirally broken pliase, obtained witli /5 < /3c ~ 1-0. Fig. ^ 
sliows tlie propagators for tt, / and PS channels on a log scale (using data obtained with 
a wall source and point sink in the latter case), resulting from approximately 40000 HMC 
trajectories of mean length 1.0. All three look to be well-approximated by straight lines, 
implying that each channel is dominated by a single particle pole. Fig. ^ shows the spectral 
functions obtained in the same three channels using Bryan's method. All three appear as 
well-localised peaks suggesting simple poles and hence stable particle states. The cross shown 
on each peak is obtained as follows. The spectral feature is fitted to a form ZG{uo — M/T) 
where G{x) is the normalised gaussian distribution, M the peak position, F the full width 
at half maximum, and Z a normalisation factor. The horizontal bar's position and width 
represent M and F respectively, and its height represents the area of ZG{uj — M/V) evaluated 
between uj—T and uj+T. The vertical error bar represents the error in this area as determined 



by the Bryan algorithm |1J]. For a narrow gaussian, of course, the central value is interpreted 
as the particle mass. 

In Table m we list the masses obtained from simulations of the U(l) model from both single 
exponential fits and MEM, as well as the area under the gaussian peak, using correlator data 
from timeslices 2-10 for the tt; for / and PS timeslices 2-8 were used. Note that for the 
lightest state, namely the tt, MEM systematically yields a lower mass, suggesting that it is 
less affected by excited state contamination, although in all cases the two methods are within 
a standard deviation. Fig. |^ demonstrates that the pion mass extracted using MEM over a 
range of bare fermion masses is consistent with the PCAC behaviour oc y/m expected for 
broken chiral symmetry. For the / and PS channels there is excellent agreement in almost all 
cases between the two methods. The PS mass is roughly twice that of the fermion, consistent 
with its being a weakly bound state. With the precision we have obtained it is possible to 
estimate the binding energy defined as Am = 2Mj — Mps; the results are tabulated in 
Table For Nf = 4 Am ^ 2.8% of the bound state mass, but the figure drops to ~ 0.15% 
for Nf = 36, which is roughly consistent with the analytical expectation that Am oc 1/A^/ 
(note, however, that the Nf = 36 simulations were performed on a smaller volume). It was 
observed in |^ that the PS wavef unction has considerably greater spatial extent for larger 
Nf, again implying it is less strongly bound. 

As discussed in Sec. |]the area under the peak is related to the strength A of the coupling 
of the operator J to the single particle state, and hence to physical decay constants. Our 
results show a systematic decrease in this coupling strength with bare fermion mass m, the 
effect being most pronounced for the vr. 

Finally in Fig. ^ we explore the effects of using different meson sources following Eqn. 
([4. 51) using data from timeslices 1-8. As in Fig. HI the spectral functions have been rescaled 



so as all to fit on the same plot. When a wall is used at either sink or source, the signal is 
completely dominated by the bound state; however, for the point-to-point correlator there 
is a significant contribution out to ua ~ 2.5. Since we have discarded data from small 
timeslices we should not expect much quantitative information from the asymptotic form of 
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p{u) in this case; indeed, as ^ oo it decays much faster than either of the ideahsed forms 
Po{uj) or pi/vi^) discussed in Sec. ^. Fig. ^ provides a graphic illustration, however, of the 
importance of choice of source in maximising the projection onto the ground state. 

5.2 Symmetric Phase 

Next we turn to the chirally symmetric phase found for (3 > f3c, where according to the dis- 
cussion of Sec. D the bound state poles should be replaced by resonances with non- vanishing 
widths. Our simulations in this section were performed for the Z2 model on a 32^ x 48 
lattice at couplings P = 0.92, 1.0 and 1.25 with 0(40000) configurations separated by HMC 
trajectories of mean length 1.0, and for U(l) on a 32^ lattice at (3 = 1.0 and 1.25 with 
respectively 30000 and 60000 trajectories of mean length 0.6. In all cases Nf = 4 fermion 
flavors were used. It proved considerably easier in this phase to simulate the model with Z2 
chiral symmetry: the U(l) simulations required a much smaller molecular dynamics timestep 
making them more expensive, and the data correspondingly of not such good quality. Data 
for the Z2 0" timeslice correlator are shown on a log scale in Fig. 0. In contrast to the broken 
phase correlators of Fig. ^ it is clear that a simple pole fit will not be successful; indeed, the 
correlators become almost flat at large t, which means that towards the centre of the lattice 
we have to worry about significant contributions from not just a backwards-propagating sig- 
nal, but also "image" sources displaced by integer mulitples of from the original source 
i- 

If we are to successfully identify spectral features as something other than simple poles, 
then it is important to study systematic effects. Fig. ^ presents results from the a channel, 
where the resonance is anticipated, showing the effects of varying the timeslice sample used 
in the MEM fit. Data from within a time window [^1,^2] were fitted; in all cases we chose a 
rather conservative value ^2 = H to minimise finite volume (actually non-zero temperature) 
effects due to the image sources discussed above, although we have checked that the results 
are insensitive to reducing t2- Fig- H shows a broad feature centred at ua ~ 0.5, whose 
"width" (actually the ratio of width to indicated by the crosses) increases as data 

from smaller times is included. Ignoring the divergence as — which we take to be an 
artifact (possibly due to a small residual vacuum expectation (cr) 7^ 0; see discussion below in 
Sec. p.3| ), the shape of the spectrum appears quahtatively similar to the large- iVj prediction 
( p.3|) . The fact that the shape of the spectrum in the massless phase is sensitive to the 
data at short times is slightly counter-intuitive, but is consistent with the observation in 
that extraction of a physical scale, namely the resonance width /i, from timeslice correlator 
data actually depends on corrections to the expected power-law falloff (yut)^^ at small values 
of fit. Note that p{uj)/uj falls to zero as ^ 00, in contrast to the constant behaviour 
expected in an asymptotically-free theory such as QCD and exemplified in Fig. |I]. The falloff 
is approximately power-law of the form a;~^, but with p fti 4 — 6, in contrast to the value 
p = 2 predicted by ( p.3|) . 

It is also legitimate to ask whether the non-zero width of the spectral feature is due 
to insufficient statistics. Fig. ^ shows the feature evolving as data is added to the sample. 
There is no significant reduction in the width of the feature as the statistics accumulate from 
0(10000) to 0(40000) configurations, although the central position and height of the peak 
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both vary slightly, supporting the conclusion that a resonance is present. 

In Fig. |l^ we compare the results from 3 different couplings. Since the artifact at u; — > 
distorts the normalisation of our result, we have rescaled each curve so that the rectangles 
of equal area to the fitted peak have the same height. The resulting curves show both 
the position and width of the resonance increasing with (3. This is consistent with ( |3.3| ), 
which predicts both are proportional to a single scale /i, if /i increases with [3 as expected. 
Within errors we find the ratio of width to central position constant and approximately 
equal to 50%. Note, however, that ignoring the spike at = the dip in the curve suggests 
lim^^o p{^)/^ ~^ 0, rather than tending to a constant as predicted by (^.31). 

Finally in Figs. and |T2| we show some results from simulations of the U(l) model. In 
this case it is possible to extract and compare spectra from both a and vr channels. The 
fitted time window is [1,10]. The bare fermion mass m is set to zero implying that for (3 > (3c 
the two channels should be physically indistinguishable, and Fig. |T2| suggests that for large 
uj this is indeed the case. There is, however, a large disparity as a; — >• between /3 = 1.00, 
where p{(jj) appears to diverge, and (3 = 1.25 where it seems to tend smoothly to zero. 
Fig. |TT] confirms that the behaviour of the correlators at large t is not really under control 
yet with the precision we have been able to obtain. Also in both cases there is more power 
in the a channel at small uj. This indicates we still lack a full understanding of systematics 
in this regime. Intriguingly, however, the large-tu behaviour is much closer to the large-A'^y 
prediction; the dash-dotted line in Fig. |12| is a fit of the form p{uj) cx uj~^'^, to be compared 
with the expected u!~^. 

To summarise, there is encouraging evidence that MEM analysis can successfully identify 
a resonance with non-zero width in this phase of the model, whose properties are consistent, 
at least in part, with theoretical expectations. Despite uncertainty about the a; — >■ limit 
that would probably require lattices considerably longer in the Euclidean time dimension to 
resolve, the MEM method is capable of yielding semi- quantitative information in this regime. 



5.3 The a Channel in the Broken Phase 

Finally we return to the chirally broken phase and switch our attention to the a channel. 
Recall that since the a is modelled via an auxiliary boson field, diagrams formed from 
disconnected fermion lines are automatically included in the calculation of the correlator. 
The main physical issues to address are whether the the a is a bound state, and if it is 
possible to detect a signal for a ^ nn decay. Conventional spectroscopy, using both simple 
pole fits and large-iVj inspired forms which include a //-threshold, have proved at best 
ambiguous for this case . Moreover due to the auxiliary nature of the field it is impossible 
to study the wavefunction, which in other channels provides clear evidence of // binding [Q. 

Fig. shows spectral functions in the a channel from simulations of the U(l) model at 
two different values of bare fermion mass m, and a comparison simulation of the massless 
Z2 model, in which there is no pion degree of freedom. We used a large statistical sample; 
respectively 1.7 x 10^ (U(l) ma = 0.01), 4 x 10^ (U(l) ma = 0.04), and 1.1 x 10^ (Z2 
m = 0) configurations were generated, and in all cases Nf = 4. Since the a has the 
same quantum numbers as the vacuum, it is necessary to subtract the vacuum term C^ac = 
X]:^,t(o"(0, 0))((t(x, t)) from the raw data to define a connected Green function. Because 



13 



of statistical fluctuations this procedure is hard to implement exactly, despite the large 
sample generated, and we believe that uncertainity in the vacuum subtraction is the origin 
of the sharp spike in the U(l) m = 0.01 spectrum centred at ua = 0.092. This feature is 
otherwise hard to explain since the lightest particle in the spectrum (see Table |I]), the vr, 
has mass rriTra ~ 0.19. We have checked that varying the subtraction constant Cyac within 
a standard deviation causes dramatic alterations to both the strength and position of this 
feature without significantly affecting the peaks at higher u, and conclude that it is not 
physical. 

Proceeding on this assumption we identify spectral features centred at ua = 0.81(2) 
(U(l) ma = 0.04), ua = 0.66(1) (U(l) ma = 0.01), and ua = 0.64(1) (Z2 ma = 0). The 
width of the features are 0(0.05) and appear stable as the number of configurations sampled 
is increased, which suggests they are not simple poles. Unhke the PS spectrum of Fig. ^ 
however, their shapes are roughly symmetric, which contrasts with the large- A^/ expectation 
that p{u) should be sharply cut off on the low-cu side but fall away more slowly on the 
high-cij side due to a // continuum. The central value of the peak for the U(l) m = 0.01 
data indicates that the state it describes is lighter than the corresponding PS state in the 
U(l) model (see Table 0), which has mass 0.77 - the // threshold in this case is at 0.793(3), 
which lies well above the point where p{uj)/u! appears to fall to zero. We deduce that for 
finite Nf there is a bound state in the a channel, which is more tightly-bound than the PS 
meson for which there are no disconnected fermion line contributions. This conclusion would 
have been difficult to reach without MEM. 

Unfortunately there is no sign of any spectral feature at the two pion threshold, expected 
following the discussion of Sec. ^ at ooa ~ 0.38 for ma = 0.01 and ua ~ 0.75 for ma = 0.04 
(implying that a — > tttt decay is certainly possible on energetic grounds in the former case). 
We have checked that there is no significant difference between spectra found in U(l) and Z2 
simulations performed at the same parameter values. Possibly this is because the height of 
the expected feature is suppressed by a power of 1/A^/ (recall (^.9])) and would need a series 
of simulations with varying Nf to expose it. Thus far, however, we are unable to report 
observation of bound state decay in this model. 

6 Summary and Outlook 

Lattice simulations of theories other than quenched QCD at zero temperature will require 
spectrum analysis techniques of greater sophistication than the currently-used single- and 
multi-exponential fits, which implicitly assume a spectral density function made up from a 
series of isolated simple poles. In this paper we have applied one of the more promising, the 
Maximum Entropy Method, for the first time to a lattice model with dynamical fermions. 
Our main findings are summarised below: 

• In the chirally broken phase of the model we have found sharply defined spectral 
features corresponding to the elementary fermion /, the simplest mesonic // bound 
state, and the Goldstone boson vr. These results corroborate earlier simulations [|], H, 
and for the first time have permitted a plausible estimate for the meson binding energy. 
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• In the chirally symmetric phase we have identified a broad resonance feature whose 
position and width agree quahtatively with the expectations of the large- A^/ approach. 
The behaviour as a; — >■ cx) is distinct from that of an asymptotically-free theory, and is 
evidence for a non-perturbative anomalous dimension associated with a UV renormal- 
isation group fixed point. 

• In the chirally broken phase we have made the first quantitative study of the a channel, 
and found that it is more tightly bound than the conventional PS meson, possibly due 
to the additional contribution of disconnected fermion line diagrams. We have been 
unable to find evidence for (j ^ tttt decay. 



Since the philosophy of the MEM method is to make the maximum possible use of data, 
we have used correlator data from as wide a time window as possible consistent with stability 
of the fit. The main problem we have faced has been systematic errors associated with the 
upper end of the time window used in the fit, particularly since we have been anxious to avoid 
finite volume effects. This has made it impossible to have control of the a; — > limit. As 
explained in Sec. |5.3| , in the a channel there may also be artifacts associated with vacuum 
subtraction. Overall, our conclusion is that MEM has proved a useful semi-quantitative 
analysis tool, but that there remains room for improvement. 

In the future it will be interesting to study 2+ld four-fermion models at non-zero tem- 
perature and/or density. Since spectral analysis requires data from many Euclidean time 
separations to be effective, it is likely to be some time before an equivalent analysis can be 
applied to QCD with dynamical fermionsQ However, in the vicinity of the deconfining/chiral 
symmetry restoring transition the dominant modification to the zero-T spectrum is expected 
to be collision-broadening due to pions, an effect absent from quenched QCD, where the light- 
est states are glueballs, but in principle present in the current model. Additionally, there is 
no longer any ambiguity about the IR cutoff, which is now T~^, and the ^ limit should 



become accessible |10]; the slope of p{uj) in this limit yields information about transport 
coefficients such as electrical conductivity via the Kubo formula |TB[ . Finally, the four-fermi 
model is currently the only model simulable at non-zero baryon chemical potential yU which 
has a Fermi surface |Q; there may be rich physics associated with phenomena such as first 
and zero sound or superfiuidity via BCS condensation to explore. 
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Figure 3: Propagators in three different channels from simulations of the U(l) model on 
322 ^ 4g ig^^^ice at (3 = 0.55, m = 0.01. 
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Figure 4: Bryan image of p{uj)/uj in three different channels using the same data as Fi^ 
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Table 1: Broken phase spectroscopy 
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Figure 5: Pion mass vs. bare mass m for (3 = 0.55, showing evidence for the Goldstone 
nature of the vr. 
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Figure 6: Bryan image of p{uj)/uj in the PS channel 32^ x 48 lattice at /? = 0.55, m = 0.01 
using correlators with different combinations of wall and point sources. 
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Table 2: Binding Energy in the PS channel 
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Figure 7: cr correlator for 3 different couplings in the chirally symmetric phase on a 32^ x 45 
lattice. 
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Figure 8: Bryan image of p{uj)/uj vs. uj in the a channel at = 1.25 on a 32^ x 48 lattice, 
showing 3 different time windows. 
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Figure 9: The same as Fig. ^ using fits from timeslices 0-11, showing the effects of varying 
the amount of data. 
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Figure 10: Rescaled Bryan image of p{uj)/uj in the a channel from timeslices - 11, for three 
different couphngs. 



23 



0.0001 



le-05 - 



le-06 - 



le-07 - 



le-08 



•--•7ip=1.00 
--7ip=1.25 

o^ap=1.00 
□--□ap=1.25 



5- 1 




-J- ^1 




^1 


- N 

- \ 


\ 




\ 

\ 

\ 


4 


\ 


/ 



-i-. 







10 



15 



Figure 11: a and tt timeslice correlators from simulations of the U(l) model on a 32^ lattice. 
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Figure 12: Bryan image of p{uS)juj in both a and vr channels using the correlator data of 
Fig. ITTl 
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Figure 13: Rescaled Bryan image of p{u!)/uj in the a channel from timeshces 1 - 10 at 
f3 = 0.65, for two different masses in the U(l) model on a 32^ x 24 lattice, and for m = in 
the Zo model on a 24^ lattice. 



25 



